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Abstract. 

The path integral approach to the quantization of one degree-of-freedom Newtonian particles is 
^^ ' considered within the discrete time-slicing approach, as in Feynman's original development. In 

the time-slicing approximation the quantum mechanical evolution will generally not have any 
stationary states. We look for conditions on the potential energy term such that the quantum 
^O ' mechanical evolution may possess stationary states without having to perform a continuum limit. 

When the stationary states are postulated to be solutions of a second-order ordinary differential 
equation (ODE) eigenvalue problem it is found that the potential is required to be a solution of a 
particular first-order ODE. Similarly, when the stationary states are postulated to be solutions of 
a second-order ordinary difference equation (OAE) eigenvalue problem the potential is required 
to be a solution of a particular first-order OAE. The classical limits (which are at times very 
^ . nontrivial) are integrable maps. 

o 

^ ' 1. Introduction 

o. 

The path integral approach to quantum mechanics was developed by R P Feynman in the middle 
of the last century |1], building on an insight of P M Dirac's j2] (see also [3]). As described 
by Feynman j3], the path integral approach "is another, a third way, of describing quantum 
mechanics, which looks quite different to that of Schrodinger and Heisenberg, but is equivalent to 
them." There have now been thousands of papers on the application of path integral techniques 
5t ! to problems in quantum mechanics. This was recently surveyed in the monograph jHj, where 
results on models that have been explicitly solved were classified and arranged in tabular form. 
In their practical application, however, path integrals are still usually performed via Feynman's 
original time-slicing approach (see, for instance, [H]). 

This Letter considers the time-sliced path integral of a non-relativistic particle with one 
degree-of-freedom, as described by a Lagrangian of the form 

Co{q,q) = le-Viq), (1.1) 

and looks for potentials, V{q), such that the quantum evolution may have stationary states whilst 
still within the discrete time-sliced approximation. 

The results have relevance to various areas of research, including the study of lattice effects 
in Monte Carlo simulations, asymptotically long-time investigations of quantum systems, and 
the investigation of discrete-time quantum systems in their own right. 
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In section 121 the potential term is derived under the assumption that the preserved quantum 
mechanical stationary states of the discrete-time evolution are eigenfunctions of a linear second- 
order ODE. 

In section El the potential term is derived under the assumption that the preserved quantum 
mechanical stationary states of the discrete-time evolution are eigenfunctions of a linear second- 
order OAE. 

2. Differential equation ansatz 

Feynman's approach gives the one-step time evolution of a wave function ip{qn,tn) within the 
time-slicing approximation as 



ip{qn+utn+i) = N exp 
Jc 



-eC{qn+i,qn) 



ij{qn,tn)dqn {2.1] 



where C is some interval of M, and n is a discrete time parameter. To obtain the wave function 
at later times the integrals are concatenated (i.e., iterated). In Feynman's original paper |3], 
the path integral is built up, and worked with throughout the paper, by employing a finite time 
interval e = tn+i — tn, which is sent to zero at the end of the calculations. In this limit (where a 
finite time interval requires an infinite number of integrals) we obtain the path integral, see |3], 
|H]. This 'time-slicing' approach is still the method that is usually employed in the evaluation 
of path integrals. As we are studying phenomena within the time-slicing approximation we set 
e = 1 without loss of generality. The discrete-time Lagrangian is (see jl] ) 

£(g„+i, g„) = -(g„+i - g„)^ - V{qn) (2.2) 

(The discrete-time Lagrangian ()2.2|) follows from the original Lagrangian (jl.ll) by the replacement 
q = (g„+i - qn)/£, see g.) For C = M, unitarity shows that |A^p = {2nhy^ (in g] iV is found 
by considering the 0th order term of an expansion in the time-step parameter). At this stage we 
refrain from stipulating the range of integration and the normalization constant N. 

A simple gauge transformation allows us to consider, instead of ()2.2j) . the (gauge equivalent) 

^Qn+l, qn) = -qnqn+1 + ^n " '^('?n) =: -gn^n+l " W{qn). (2.3) 

In this section, we look for forms of the potential, V{q), (or, equivalently, W{q)) such that 
after evolution by equation (|2.ip the new wave function is a solution to the same eigenvalue 
problem of a second-order ODE, with the same eigenvalue (and so will be for any number of 
evolutions). We write the eigenvalue problem as 

M,^{q) = E^{q) (2.4) 

where Mq is an arbitrary linear second-order ODE (which for convenience is written as a partial 
differential operator) 

Mg=Po{q)g^+Pi{q)-g-+P2{q). (2.5) 
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The new wave function, evolved by (|2.H) . although required to be a solution of the same eigenvalue 
problem, is not necessarily the same function. However, the requirement of unitarity will remove 
the possibility that the evolved wave function becomes the identically zero solution. 

Over-tildes will now be used to denote time updates, i.e., q = g„_i, q '■= Qn, Q '■= In+i, etc. 
It is required that 



However, by definition 



(2.6) 



E^{q) 



N I exp 

'c 



^>C(g,g) 



c 



where S is the surface term 



A^ / <^Mgexp 



^>C(g,g) 



Eil){q)dq 

ip{q)dq 



N I exp 

c 



^>C(g,g) 



MMq) \dq 



S, 



S 



Po(g)exp 



h 



C-iq^q) 



^#(g) + -^) 



(2.7) 



(2i 



c 



and Mq is the adjoint of Mq. For ()2.7|1 to hold it is required that the integral in equation p.lj) 
exists, and, furthermore, the integral converges uniformly with respect to q in an appropriate 
region and the integrand of the integral relation is an analytic function of q in this region when 
q is on the line of integration. It will be assumed that these analytic requirements hold. The 
surface term vanishing, i.e., S" = 0, seems to be generally in accord with physical requirements 
(for instance, if C = M, then S* = is achieved by requiring that V'(?) ^"^^ 
\q\ — * oo). If the surface term vanishes then it is required that 



-|^ tend to zero as 

dq 



M„ 



exp 



h 



C-iq^q) 



Mq exp 



n 



>C(g,g) 



(2.9) 



This PDE is a functional equation for both Mq and £(g, q). It can be solved as follows. 

Substituting the form of the Lagrangian ()2.Hj) into ()2.9j) and comparing coefficients of 
different powers of q reveals Po{q), Pi{q) and P2{q) are quadratic polynomials. Hence, set 



-^Po{q) = Aof + Boq + Co 



-tPM =Aif + Biq + Ci 

h 



P2iq) = A2r + B2q + C2 



and substitute this back into ()2.9|) to see how many of the coefficients are actually independent. 
A comparison of the g^ terms of ()2.9|) shows 



Aoq^ + Aiq + A2 = A^q^ + B^q + Cq. 

Comparing coefficients of different powers of q reveals Ai 
A comparison of the q terms of ()2.9p shows that 

Dq^ + Eq + F 



(2.10) 



Bq and A2 = Cq. 



W'{q) 



Aoq^ + Boq + Cq 



(2.11) 
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where we have introduced 

D:=^{Ai + Bo) E:=B,-th2Ao F ■=^{B2 + Ci) - ihB^. (2.12) 

A comparison of the (f terms of (j2.9|) is then made: The g° terms reveal that for (j2.9|) to 
hold, either B2 = Ci oi F = Q. As F = limits the form of the potential term, B2 = Ci is 
chosen. The q^ terms reveal that for ()2.9|1 to hold D = Ai. However, Ai = Bq and therefore 
D = Bq. The q^ terms are trivial. 

Taking the above findings together gives the form of Mq. For aesthetic reasons we make the 
replacements a := Aq, (3 := i?o, 7 := C'o, e := E, ( := F, and set C2 = ihe/2 so that Mq is 
formally Hermitian, to obtain 

Mq = ~h\aq^ + /5g + 7)77^ + iKPQ^ + (e + 2zna)g + ihp + C)^ + W + (^^/^ + Qq + ^-(2.13) 

ag^ oq I 

The potential term follows from the expression for Vr(g), and we have 

^'(^) = ^T^^^ where a,/5,7, e, C e M. (2.14) 

otq^ + pq + 'y 

The quantum evolution given by equations ()2.1|) . ()2.3|) and ()2.14|) (which is defined up to an 
unimportant constant of unit magnitude) sends the eigenfunctions of (j2.4p . where Mq is given by 
()2.13|) . to eigenfunctions of the same eigenvalue problem ()2.6|) with the same eigenvalue, if the 
analytic requirements given between equations ()2.7|) and ()2.9|) hold. If the singularities of the 
rational function W'{q), (|2.14j) . (which are the singularities of the ODE given by (j2.13p ) occur 
on the real- line then there are also issues to be addressed in the integration of W\q) to obtain 
W{q). (These are issues associated with the subset of M on which the model is defined, which 
may have singularities of W'{q) as boundaries. There are numerous possibilities here, which must 
be dealt with on a case by case basis.) 

With the parameters in generic position, the finite singularities of the ODE given by ()2.13p 
are regular. According to the classification scheme of jT^j they have s-rank= 1, while the 
singularity at infinity is irregular with s-rank=2. Therefore the ODE given by (j2.13|) has s-rank 
multisymbol {1, 1, 2}; hence it is of confluent Heun class. As a consequence, the analytical work 
required to show that the quantum evolution sends an eigenfunction to the same eigenfunction 
will be quite involved in general. One subcase has already appeared in the literature; namely 
the case where a = /3 = C = and C = R. This is the discrete-time harmonic oscillator, p]. In 
exact analogy with the continuous-time quantum harmonic oscillator, the eigenfunctions are of 
the form Hk{aq)e~^ ^ , a, 6 G M, where if^ denotes the fcth Hermite polynomial. All analytical 
requirements are met in this case, and the path integral evolves eigenfunctions to the same 
eigenfunctions. That is, these are quantum mechanical stationary states of the system whilst 
still within the time-slicing approximation. 

The classical limit of the quantum system given by equations ()2.1|) . ()2.3p and ()2.14|) is 
an integrable map of the plane known as the McMillan map J2] (see also |7j). The classical 
limit of the invariant is obtained using —ih-^ = p. In this connection it should be noted 
that integrable quantum mappings, in terms of Heisenberg operators, possessing exact quantum 
operator invariants have been proposed before, e.g, in [T3], [Ej. 
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3. Difference equation ansatz 

We now consider the quantum mechanical evolution associated with forms of the potential, V{q), 
(or, equivalently, W{q)) such that the wave function is a solution to an eigenvalue problem of a 
second-order O AE both before and after its time-evolution within the time-slicing approximation. 
Furthermore, the original and evolved wave functions are associated with the same eigenvalue. 
We write the arbitrary linear second-order OAE eigenvalue problem as 



Nqi^iq) = Ei){q) Ng = pi{q) exp 



d_ 

dq 



+ Po{q) +p_i(g)exp 



d_ 
dq 



Define, for convenience, w{q) := e n^^'i\ It is required that 



Ngi^iq) = E^{q). 




, by definition 


Etjj{q) = N exp 
Jc 


{m,) 



Eip{q)dq. 



(3.i: 



(3.2) 



(3.3) 



It will be assumed that the integral in ()3.3|) exists, and, furthermore, that it converges uniformly 
with respect to q for g G M. It is also assumed that the integrand is an analytic function of q for 
all g G M. In this section the usual physical assumptions C = R and iplq) G L^(M) are made, so 



oo 

exp 

-oo 


i ^ 


oo 

exp 

'OO 


i 



EiP{q) = N 

= N 

+ w{q + a)p_i(g + a) exp 
Hence, for fj3.2|) to hold, it is required that 



w{q)\pi{q)'ilj{q + a) + Po{q)^/j{q) + p^i{q)ilj{q - a) \dq 



< w{q — a)pi{q — a) exp 

ilj{q)dq 



r' 



+ w{q)po{q) 



-r' 



(3.4) 



w{q — a)pi{q — a) exp 



r' 



+ w{q)pQ{q) + w{q + a)p_i(g + a) exp 



-r'' 



w{q) exp 



-r' 



pi{q) + w{q)po{q) + w{q) exp 



r"^ 



p-i{q). 



(3.5) 



This is a functional equation that is solved in a similar manner to the development of the previous 
section to obtain pi{q), Po{q), P-i{q) and w{q). (In this case using the independence of different 
powers of e^"'^, rather than q.) 

After some tedious but straightforward calculations the following first-order analytic OAE 
for W{q) is obtained. 



W{q + a)-W{q) = ih\n 



Ai exp [|ag] + Bi + Ci exp [-fag] 
A_i exp [^aq + {a"^] + Aq + Ai exp [-{aq - j:a'^] 



.(3.6) 
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The difference operator (|3.1|) . of wliicli tlie postulated quantum meclianical stationary states are 
eigenfunctions, is 



N. 



Aiexp 



h 



aq 



+ Bi + Ci exp 



h 



aq 



exp 



d_ 
dq 



+ Aq exp 



+ < A_i exp 



h 



aq 



+ Ao + Ai exp 



h 



aq 



exp 



— a- 



d_ 
dq 



aq 



+ Bo + Bi exp 



h 



aq 



(3.7) 



Tlie aim of this investigation is to find potentials of one degree-of-freedom Newtonian 
systems such that their quantum mechanical time evolution has stationary states even within 
the time-slicing approximation. Hence we wish to find explicitly real-valued expressions for 
W{q) (and, consequently, V{q)). To this end, the following parameter replacements are 
made, A_i = ie-*4(5 + /? + 2ia), Aq = i(C + 27), ^1 = jer.4(S - /5), Si = i(C - ^7), 
Ci = ^e~"h~(6 + P — 2ia); and substituting into (j3.6p gives 

'/3sin[|(g + f)]+acos[|(g + f)]+7^ 



W{q + a) -W{q) = 2htan 



{3.1 



.5cos[t(g + f)]-asin[|(g + f)]+C. 
The difference operator Nq is also easily rewritten in terms of the new parameters and 
trigonometric functions. 

The first-order equation obtained by exponentiating ()3.6|) is easily solved, in a formal 
manner, by a product of four hyperbolic gamma functions ^^, fZj. Ruijsenaars' hyperbolic 
gamma functions \1^ are so-called 'minimal solutions' that are essentially unique, and have 
distinctive asymptotic properties (the w{q), of course, is only defined up to a factor with period 
a) . Alternatively, in the difference case one may consider a scenario where the path integral is set 
up in terms of discrete summations, or Jackson g-integrals ^Ul , and where the first-order O AE, 
(I3.fi|l or ()3.8|1 . is solved by an iterative procedure. In a specific case there may be grounds for 
choosing a particular solution approach. This, however, falls outside of the remit of the current 
Letter, which stops at giving the equation that the potential must obey (i.e., ()3.6|) or ()3.8|) ). 

The classical limit of this system is found after first making the reverse Schrodinger 
replacement —ifi^ = P (and p = —q in this gauge). However, in a naive classical limit Nq 
trivializes. On making the replacement a = ash, where as is independent of h, the system now 
has a rich classical limit, namely the (integrable) Suris map ^Hl- The passage to the classical 
limit is nontrivial; the equation that the potential must obey is qualitatively changed by the 
process: the first-order difference equations that define W{q) (either (j3.6|) or, equivalently, (j3.8|) ) 
become first-order differential equations. From an alternative perspective to that of this Letter, 
it can be remarked that we have an 'integrable quantization' of the Suris map, where, to retain 
integrability, the potential term has nontrivial quantum, 0{h), corrections. (In other settings it 
has long been realised that to retain integrability after quantization the potential term may need 
quantum deformations, see [Hj.) 

It is satisfying to note that in both cases the functional equations f ()2.9p and ()3.5|) ) give, in a 
straightforward manner, the specific form of the second-order differential or difference operator, 
as well as the equation for the potential term, such that the discrete-time quantum evolution 
may possess stationary states. 
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